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Abstract 

The dynamics of a tachyon field plus a barotropic fluid is investigated in spatially curved FRW 
universe. We perform a phase-plane analysis and obtain scaling solutions accompanying with a 
discussion on their stability. Furthermore, we construct the form of scalar potential which may 
give rise to stable solutions for spatially open and closed universe separately. 
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I. INTRODUCTION 



The nature and the origin of dark energy is a fundamental puzzle in modern cosmology. 
Most dark energy models can be constructed by using a slowly rolling canonical scalar field, 
termed as quintessence. However, there has been increasing interests in alternative models 
with a non-canonical kinetic term. Among these models the most general formalism perhaps 
is k-essence A more specific case is the tachyon field 2I, which is motivated from string 
theory. It can be viewed as a special case of k-essence with Dirac-Born-Infeld (DBI) action 
Although the tachyon is an unstable field, its state parameter in the equation of state 
varies smoothly between —1 and 0, thus many authors have already considered the tachyon 
field as a suitable candidate for a viable model of dark energy phenomenologically j4-13|. 
For a review, we can refer to Ref . fl^ . 

However, the dynamical dark energy models driven by a scalar field suffer from the so- 
called fine-tuning problem and coincidence problem. In order to address these problems, one 
may employ scalar field models exhibiting scaling solutions , [iSHSJ] . The scaling solu- 
tions as dynamical attractors can considerably alleviate these two problems. Furthermore, 
by investigating the nature of scaling solutions, one can determine whether such behavior 
is stable or just a transient feature and explore the asymptotic behavior of the scalar field 
potential. 

Many authors have investigated a lot of scalar field models containing scaling solutions. 
For instance, a canonical scalar field with an exponential potential has scaling attractor 



solutions 



16] . For quintessence dark energy model, there are two scaling solutions. One 



is fluid-scalar field scaling solution, which remains subdominant for most of the cosmic 
evolution. It is necessary that the scalar field mimics the background energy density (radi- 
ation/matter) in order to respect the nucleosynthesis constraint and can also alleviate the 
fine-tuning problem of initial conditions. The other is scalar field dominated scaling solu- 
tion, which is a late time attractor and gives rise to the accelerated expansion. Since the 
fluid-scalar field scaling solution is non-accelerating, we need an additional mechanism exit 
from the scaling regime so as to enter the scalar field dominated scaling solution at late 
times. For the discussion on the exiting mechanism, we can refer to Refs. 



For tachyon field dark energy, the scaling solutions have also been investigated by many 
authors, for example Refs. {4- 13 1. To be considered as a realistic model of dark energy, it is 
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found that the fluid-scalar fleld scahnff solutions are absent and only the scalar field domi- 

n 

nated scaling solutions exist [6|. This is very different from the quintessence case. Therefore, 
just as pointed out in 13|, tachyon models require more fine-tuning to be consistent with 
observations. Nevertheless, it is worthwhile to point out that the fiuid-scalar field scaling 
solutions may be obtained when the Gauss-Bonnet coupling between tachyon field and fiuid 
is considered in the tachyon dark energy model 34 1. 

Although the latest results of WMAP5 have placed a constraint, —0.063 < Qk < 0.017 



40| . on the fiatness of our observable universe, indicating that our observable universe is very 



close to fiatness, it is still possible that our observable universe is spatially curved. Therefore, 
it is also interesting to investigate the dynamical behavior of dark energy models in a spatially 
curved FRW universe. Recently Copeland et.al. have extended such investigations to the 
quintessence model in spatially curved FRW universe 4l)] , following the strategy they have 
developed in Ref. 28 1. Sen and Devi [ll| have also explored the scaling solutions with tachyon 



in modified gravity model employing the same method. In this Letter, we will closely follow 
this route to investigate the dynamics of tachyon dark energy model in spatially curved 
universe. 

Our Letter is organized as follows. In Section 11 we present the associated equations of 
motion for the tachyon field including the background fiuid and obtain the scaling solutions. 
Then, we analyze the stability of these solutions. In Section III we turn to construct the 
scalar potential leading to such scaling solutions. In particular, its asymptotical forms are 
obtained in various circumstances for spatially open and closed universe respectively. 



II. DYNAMICS OF THE TACHYON FIELD IN THE PRESENCE OF A FLUID 

Let us start with a cosmological model in which the universe is filled with a tachyon field 
evolving with a positive potential V{(j)) and a barotropic fiuid with an equation of state 

= (7 — l)p^, where 7 is the adiabatic index. We note that 7 = 1 for a pressureless dust 
and 7 = 4/3 for radiation. The pressure and the energy densities of the tachyon field are 
respectively given by 

P4> = -v{<p)^/l-^^ (1) 




(2) 
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then the effective adiabatic index of the tachyon field is given by 



1^ = = • (3) 



From the above equation, we can see that < 7^ < 1 



As shown in Ref . [28|(4l| , either in a spatially curved FRW universe or a flat universe 
but described by modified gravity models, the Friedmann equation of the universe can be 
uniformly written as an effective form 

= P-L\p{a)), (4) 

where we have set 'i'K jm^ = 1. H = a/ais the Hubble parameter and p(a) is the total energy 
density of the universe. In this paper, we adopt the notion H = eV H"^, where e = +1 for 
expanding universe and e = —1 for contracting universe. The function L{p{a)) is assumed 
to be positive-definite without loss of generality. We also note that when L{p{a)) = 1, the 
effective Friedmann equation (j4]) can be reduced to the standard spatially fiat case. 
The energy conserved equation of matter and tachyon field are respectively 

+ S'y Hp^ = 0, (5) 
^ + 3H^+^=0. (6) 



1 - (^2 - V 

Now we intend to demonstrate a phase-plane analysis on this system. As usual, we define 

X^^, Y.^, (7) 
VP 

where p is the total energy density of the universe and the subscript N denotes the derivative 
with respect to the number of e-folds, = Ina. As a result, Eq.(j4]),(|5])and(l6]) can be 
rewritten in the form 

X^N = {X^ - 1)(3A: - eVSYX), (9) 
y„ = jl-eVTiXYX + 37 - —^^=^], (10) 
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and the constraint equation for the total energy density becomes 

, + ^ = 1. (12) 

Since < -^^^ < 1, the allowed range of X and F is < + < 1, namely < < i 
and < X2 < 1. 

We consider the special case A = const. Taking a derivative with respect to the scalar 
field 0, then from Eq.(IH]) we obtain 

„ 3 dlnL 



The fixed points for this system can be obtained by setting X = and F^v = in Eqs. 
dH]) and (fTOj) . Essentially we have four fixed points: 

X, = 0,F, = 0, (14) 

Xe = ±l,F, = 0, (15) 

X, = ^Ys,Y, = Ys, (16) 



where F^ is defined by 



X, = ±V7,F, = ±^, (17) 



We now investigate the stability around the critical points by evaluatin g th e eigenvalues 
of the matrix Ai. The way of evaluating eigenvalues has been given in jl^|l6| 42 1. Here we 
list our analysis corresponding to each solution as follows and then summarize the results 
in Table I. 

1. For the fluid dominated solution(Xc = 0, Fc = 0), the eigenvalues are pi = —3,fi2 = 
37/2. Therefore, this critical point is an unstable saddle point for 7 > 0, whereas it 
is a stable node for 7 = in an expanding universe (e = 1). However, since the fixed 



point leads to the density parameter fi,^ = ^ = 0, it cannot be used as a late-time 
attractor. 



2. The fixed points {Xc = ±1 , Fc = 0) correspond to the scalar field kinetic dominated 
solution with 7^ = 0^. The eigenvalues are fii = 6, fi2 = 37/2, indicating that 
this point is unstable in an expanding universe. However, in a contracting universe 
(e = —1), it implies unconditional stability when these solutions exist. 



3. For the well-known solution (X^. = XYs/\/3, Yc = Yg), we have eigenvalues fii = 
—3 + ^Y^, fi2 = —37 + A^F/, where satisfies — 3 < /ii < —3/2. Moreover, we have 
/i2 < for 7 > 7s = ^^s^ indicating that, in an expanding universe, the fixed point is 
stable for 7 > 7^, whereas it is an unstable saddle point for 7 < 7^. On the contrary, in 
a contracting universe the critical point is unstable for 7 > 7s and be a saddle point for 
7 < 7s. Since in this solution, the kinetic energy of the tachyon field is proportional to 
the potential energy, we call this solution as the scalar field dominated scaling solution 
as Ref.j4l|. 



4. The value 7^ at the last critical points (Xc = ±^7, Yc = is 7</, = 7, which 

means both energy densities and Py decrease with the same rate. The eigenvalues 



are /ii,2 = f [7 - 2 ± ^1772 - 2O7 + 4 + §-f^^/T^^. From Eq.(^^ we have < 7 < 
7s = ^Yj^, note that 7s is always smaller than 1. The real parts of /ii and /i2 are 
both negative under the condition 7 < 7s. Obviously, it is an unstable solution in a 
contracting universe, whereas in an expanding universe the scaling solution is always 
stable. However, we need to caution that the existence of the scaling solution requires 
the condition < 7 < 7s < 1, which is not satisfied for known realistic fiuids. So it is 
always treated as an unpractical solution. 
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A > 
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e = -1 


VA, V7 
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exists 


VA, V7 


when 
exists 


VA, V7 


No 


A > 

< 7 < 7s 


No 


A < 

< 7 < 7s 


No 



TABLE I: The existence and stability conditions for an expanding (e = 1) and a contracting 
(e = — 1) universe, containing a tachyon field and a fluid with the adiabatic index 7. 



In addition, we also plot some representative figures (FIG.l and FIG. 2) for the evolution 
of universe, from which we can see the better picture of the late time acceleration. From 
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FIG. 1: The evolution of the parameters x and y in the presence of a barotropic fluid with 7 = 1. 
Here we take A = 1 and the initial conditions Xi = 0.8, Ui = 5.0 x 10"'* when = 0. 

FIG.l, we can see that when 7 = 1 > 7^, the evolution of the parameters x and y approach 
the critical point {X^ = AY,/a/3 ~ 0.53, = ~ 0.92) as N increases. From FIG.2, 
we can see that in the attractor regime, the equation of state ^ —0.72 and the density 
parameter Q^n = but Qfp = 1 which is just tachyon field dominated epoch. 

III. TACHYON FIELD POTENTIALS IN CURVED UNIVERSE 

In the rest of our Letter we focus on the scaling solution given by Eq.( fT6l) . and construct 
the tachyon field potential which may give rise to this late time attractor, provided that 
the specific forms of L{p) is given for spatially curved universe. At critical points we always 
have Yc = ^ = constant, thus L{p) can be described as L{V). Substituting this into Eq. 
([8]), we can obtain 



Therefore, given a specific form of L, it is possible to derive the potential for tachyon 
field in the scaling solution by integrating Eq.( IT9l) . 




(19) 



7 




N 



FIG. 2: The evolution of 0,^, ilm, w^p, and for 7 = 1. Here we take A = 1 and the initial 
conditions Xi = 0.0085, yi = 0.0085 when iV = 0. 

A. Open FRW universe 

In open FRW universe, L(p) is given by 



L{p) = Jl + —,, (20) 
V pa^ 

At first we assume A 7^ and discuss the form of the potential in an expanding universe. 
Then, the special case of A = is presented at the end of this subsection. 
■ case (i) A ^ 

When the universe is expanding (e = 1), the solution given by Eq.( !T6|) is stable for 
7 > 7j, = |A^y^^. By using the energy conserved equation of tachyon field 

P0 + 3-7<^p0 = 0. (21) 
a 

Therefore, we can obtain 

« = a«Pr^P;'^- (22) 
Because this solution is the scalar field dominated scaling solution, the fiuid is absent. 
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Therefore, we can have 



a = a^,)p'''(^)p-^. (23) 
With the help of the expression of the effective adiabatic index = = ^X^Y^, the 



scale factor a is given by 



1 



a = Ap ^,A = a(,)P('.)''= , (24) 

where the subscript (i) refers to some initial time. 

The correction function L given by Eq.( l20l) is then rewritten as 



Lip) = \l^ + j^P^, (25) 

where p = — 1, and for the region of existence of these solutions, the valid range of fj, 
is —1/3 < p < oo. 

At the critical point, is a constant. Substituting Ys = \/V/p into Eq. (fT9|) . the scaling 
solution potential will correspond to 

dV 

= -A0. (26) 



1 + jrYs-^^Vt^ 



Subsequently, we will classify the asymptotic behavior of the scalar field potential above 
in terms of the sign of the parameter /i. In addition, we also note that the equation f l26p 
is invariant under the transformation A — )■ —A and (f) — )■ —(p. Therefore, our work can be 
restricted to the first quadrant without loss of generality. 

1. For fi < [i.e.X'^Yg > 2), in an expanding universe at early times, where the curvature 
is negligible, the potential for the tachyon field is obtained a.s V ^ 4A~^0~^. However, 
at late times, the universe becomes one dominated by the curvature, the asymptotic 

2 2 

form of the potential is a function like oc A m+i . 



2. For yU > {i.e.X^Y^ < 2), we find that at early times, the potential relating to the 

2_ 2_ 

curvature dominated universe can be approximated by oc A a^+i m+i . Once the 
curvature is negligible at late times, it can be shown that the potential is ^ ~ 4A~^0~^. 

3. For the special case of p = 0, the correction function L becomes a constant, and 
it can be seen as a modified Newton's gravitational constant in Eq.dl]). The evolving 
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background corresponds to a flat Friedmann universe, and the potential asymptotically 
has a form of the late-time attractor V ~ ^^^0"^. Note that in this special case, the 
scalar fleld and the curvature scale together. 

We summarize above results in the Case A category of Table II. 

However, according to the results of Table I, we know that in a contracting universe, the 
scalar fleld dominated solution is not stable, but the kinetic dominated solution is stable. 
We must point out that in kinetic dominated solution, since = 0, we have the potential 
Vi<P) = 0. 

■ case (ii) A ~ 

For the special case of A = 0, the potential will be a constant. The solution described 
by Eq.( |T6l) is still applicable and from it we obtain Xc = 0, Yc = 1, which is nothing but 
a de-sitter solution, as substituting this into Eqs.([T])and ([2]) we have = and the state 
parameter a; = — 1. The exact solution can also be obtained by substituting these into the 
Eqs.dH), © and ©. 

From Table I, we can conclude that the tachyon fleld dominated solution is stable in the 
expanding universe, but unstable in a contracting universe. 

B. Closed FRW universe 

Now we turn to obtain the tachyon potential at the stable critical point for a closed FRW 
universe. In this case L{p{a)) is given by 

L{p) = Jl-^- (27) 
V pa 

Obviously, the valid range of the correction function is 1 > L > 0. Specially when 
L = 0, the scale factor is a constant, which neither expanding nor contracting. Therefore, 
we can assume at that point, the universe probably experiences a bounce and the evolution 
transmits from expanding(contracting) to contracting (expanding). Let us follow this idea 
to do the rest of the analysis. As before, we begin our discussion on the expanding universe 
and then analyze the difference between expanding and contracting universe for the scaling 
solution. 

▼ case (i) A 7^ 
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In an expanding universe(e = 1), the scaling solution given by Eq. (fT6|) exists for 7 > 
7s = —3-^, and the scale factor is still given by Eq.f l24p . Then we can rewrite the correction 
function as 

Lip) = \J 1 - -^P'- (28) 
Subsequently, using Yc = ^ and substituting (l28l) into Eq. (fT9l) . we obtain 

/ , , = -A^. (29) 

As the case in the open FRW universe, we discuss the asymptotic behavior of the scaling 
potential by classifying the valid range of /i into positive, negative and vanishing regions. 
We can also choose to work in the first quadrant without loss of generality. These are 
summarized in Case B of Table II. 

1. For fi > (i.e. A^F/ < 2), at early time of an expanding universe, we have 
a'^^'^s < ^-2^ However, the universe cannot be dominated by the curvature as the 
correction function is bounded in a closed universe, otherwise the right-hand side of 
the Priedmann equation @ would be negative. Consider a pragmatic point, the curva- 
ture is just subdominant to the scalar field. When L(p) — )■ 0, the total energy density 
has a maximum p — )■ Pmax = i^Y^^^ and the universe has a constant scale factor. 
From Eq.(I71), we conclude that the potential is almost a constant with V ~ Y^(^y^^. 
Subsequently, the universe starts to expand, the tachyon field gradually become dom- 
inated over the curvature. At last, one can neglect the contribution of the curvature. 
Therefore, the universe can be considered a fiat FRW spacetime, and the asymptotic 
form of the potential become V ~ 4A^^0^^. 

2. For yU < (i.e. X'^Y^ > 2), the universe dominated by the tachyon field starts to expand 
from an approximate flat FRW spacetime, and the scaling solution potential has the 
form V ~ 4A~^0~^. As time goes on, the total energy density will reach the maximal 
value, p — )■ Pmax = i^y^'^, where the contribution of the curvature and the scalar 
field are almost equal. The scale factor does not change as -ff — )■ 0, and the potential 
is almost a constant V ~ Y^{^y^'^. After this turning point, the universe changes 
from expanding to collapsing. From Table I, we see that in a contracting universe, the 
late time scalar field scaling solutions are unconditionally unstable. Once the universe 
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starts to collapse, the solution will asymptote towards the kinetic dominated solution 
with 7<^ = 1. 

3. For /i = 0, the situation is similar to the open universe scenario. The correction 
function L becomes a constant. The potential has the form V ~ where 
the contributions of the scalar field and curvature to the Friedmann equation scale 
together. 

As mentioned in the subsection for open universe, for a contracting universe if the kinetic 
dominated solution exists they will be unconditionally stable, and the potential V{(f)) = 0. 
▼ case (ii) A ^ 

Similar to our description of an open universe scenario, this case corresponds to a constant 
potential, and only the fixed point described by f|T6l) exists. For A = 0, which results in an 
expanding universe, the attractor solution reduces to = and = 1, describing an exact 
de-Sitter solution. We notice that the assumption of taking as a monotonically varying 
function of time breaks down if oc = 0. 

Furthermore, based on our previous analysis it is worth noting that when A ~ 0, the 
critical point f lT6|) is unstable in a contracting universe, but becomes a stable one in an open 
universe. 

Finally, from the above analysis, we note that the evolution of the universe including 
the curvature term is same as the case without the curvature term 6|]. However, due to 
the curvature term, the required potential is different. As shown in {g], in the case without 
curvature term, the required potential is always V{(j)) ~ 0~^. But when the curvature term 
is taken into account, the potential has various form in the different epochs of the universe. 



IV. CONCLUSION AND DISCUSSION 

In this Letter we have investigated the dynamics of a tachyon field in FRW universe with 



spatial curvature. Following the scheme presented in [28|, |4l| , we denoted the modification of 
the cosmological equation due to the spatial curvature by a general function L{p), and then 
derived the conditions under which the system can enter a scaling solution. In particular, we 
obtain an attractor solution, where the tachyon field dominates over the fluid and the kinetic 
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Case A 


Case B 






;U = 


/i > 


< 


;[i = 


/X > 


Early times 


V ~ 

negligible 
curvature 


V ~ 
scaling 
curvature 


2 

curvature 
dominated 


negligible 
curvature 


V ~ ct)-"^ 
scaling 
curvature 


curvature 
only just 
subdominant 


Late times 


2 

V 4> f +1 

curvature 
dominated 


V ~ 
scaling 
curvature 


negligible 
curvature 


curvature 
only just 
sub dominant 


V ~ 

scaling 
curvature 


V ~ 

negligible 
curvature 



TABLE II: This table summarizes the asymptotic behavior of an expanding (e = 1) universe 
described by the scaling solutions. Case A refers to open FRW universe, and Case B corresponds 
to closed FRW universe. 



energy of the field scales with its potential energy. Furthermore, given the modification 
function for an open and closed universe respectively, we discussed the form of the scalar 
potential which gives rise to the late time attractor solution. For spatially open universe, 
we conclude that in regions where the curvature is negligible, the asymptotic form of the 
potential is oc 0"^^ while in regions where it is dominant, the approximated potential will 

2 

heV (X (j) f'+i . As far as the spatially closed universe is concerned, we find that the curvature 
is not allowed to be dominant among the total energy density of the universe, but only 
subdominant to the scalar field density, where the potential is a constant as \^ ~ Y^^^Y^^ . 
When the curvature becomes negligible, the potential can be approximated by ~ 0"^ 
Moreover, it was noticed that in a contracting closed universe, the scalar dominated solutions 
are not stable, and the kinetic dominated solution will be a late time attractor if it exists. 
Comparing with the quintessence scalar field dark energy model with spatial curvature 



4l| . the fluid-scalar field scaling solutions is absent in the tachyon dark energy model. In 
addition, in the quintessence scalar field dark energy model, there are the two forms of the 
potential: the exponential potential and the power-law potential. But in the tachyon dark 
energy model, only the power-law potential is required. 

In this Letter we only consider the case that the scalar field and the fluid do not interact 
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with each other. It is a very interesting question if such an interaction term can be introduced 
in order to obtain a fluid-tachyon field scahng solutions in curved FRW universe. Effort has 
been made in present framework, unfortunately, the answer seems to be negative. However, 
it should be very desirable if we consider adding a Gauss-Bonnet coupling as in {34]. We 
expect to make further progress along this direction. 
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